August 18, 2008 



The Dynamics of the One-Dimensional 
Delta-Function Bose Gas 

Craig A. Tracy 

Department of Mathematics 
University of California 
Davis, CA 95616, USA 
email:tracy(Smath.ucdavis . edu 

Harold Widom 

Department of Mathematics 

University of California 
Santa Cruz, CA 95064, USA 
email: widom(Sucsc. edu 

Abstract 

We give a method to solve the time-dependent Schrodinger equation for a system 
of one-dimensional bosons interacting via a repulsive delta function potential. The 
method uses the ideas of Bethe Ansatz but does not use the spectral theory of the 
associated Hamiltonian. 
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1 Introduction 



In this paper we give an alternative approach to solve the time- dependent Schrodinger 
equation for a system of one- dimensional bosons interacting via a delta function potential. 
This quantum many-body model, called the Lieb-Liniger model [6l[7], is "exactly solvable" 
by means of Bethe Ansatz [2] in the sense that the ground state energy |7] , the excitation 
spectrum [H| as well as equilibrium thermodynamics [13] are known. (See [5], [TU] for 
textbook treatments.) Current widespread interest in the Lieb-Liniger model has arisen 
because of its connection to ultracold gases confined in a quasi one-dimensional trap. 
Indeed, it has been recently shown that the Lieb-Liniger model for one-dimensional bosons 
"can be rigorously derived via a scaling limit from a dilute three-dimensional Bose gas with 
arbitrary repulsive interaction potential of finite scattering length" [S] . These applications 
focus interest on the time-dependent solutions. For a review of these developments see 

Recall that the Lieb-Liniger 5-function (A^-particle) Bose gas Hamiltonian is 

j=l ^ j<k 

where c > (repulsive case) and xj G mQ Most work P H El [71 El [12] focuses on the 
eigenfunctions and eigenvalues of H from which, in principle, any solution \l/(x;t) of the 
time-dependent Schrodinger equation 

(2) 

subject to the initial condition 

^(x; 0) = ^(xi, . . . , Xiv; 0) = ^q{xi, . . . , xjv), (3) 

may be constructed!! In practice, due to the complexity of the Bethe eigenfunctions and 
the nonlinear Bethe equations, it is difficult to analyze time-dependent solutions whose 
initial conditions are not eigenfunctions. We give here an alternative to the spectral 
method for solving ([2]). 

Since different physical conditions require different choices of the initial wave func- 
tion t/'o, we develop here a flexible mathematical structure to incorporate these different 



-'^We choose units where 2m = 1 and h = 1. 

^V'o is a given symmetric function of the coordinates Xj. 
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choices. Suppose we solve ([2D with the initial condition 



1 ^ 

M^) = M^i^ • • • ' = ]^ 5Z n (^'^o) ~ y^) 

o-eSjv j=i 

where Sat is the permutation group acting on {1, . . . ,N}. Here yj G M are fixed and 
without loss of generality we may assume 

yi < y2 < ■ ■ ■ < yN- (5) 

If \E'5(x,?/;t) denotes this solution, then the solution to (E]) satisfying ([3]) is 

= / '^six,y;t)My)dy 



{dy = dyi ■ ■ ■ dy^). The subject of this paper is a Green's function for ([21). 

We remark that there are two natural choices for the initial wave function ipQ. The 
free particle Gaussian wave function localized at y i^ 



Thus an initial condition where the particles are separated and free is 

N 

^free(a^) = c^v ^1 Yl'fa{Xaijy,yj,Pj)- 

where cat is a normalization constant. 

A second choice is for the Bose gas to be initially confined to a subset of M and in its 
ground state. Solving ([2]) with this initial condition corresponds to the confined Bose gas 
freely expanding into all of M. The ground state wave function for confinement to M"*" has 
been computed by Gaudin ^ and confinement to a hard wall box by Batchelor, et al. [1] . 

We now recall a well-known [7], [10] reformulation of the problem. Since we seek 
symmetric solutions, it is sufficient to solve ([2]) in the region 

IZ : —00 < Xi < X2 < • • ■ < Xiv < 00. (6) 



■^Thus \ipa{x;y^p)\^ is the Gaussian distribution centered at y with variance a?. Solving the time- 
dependent free Schrodinger equation shows that {x) = 2pt\ and hence, 2p is the classical velocity. (Recall 
2m = 1.) 
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In the interior of TZ, e.g. Tl° : Xi < X2 < ■ ■ ■ < xn, the (5-functions are zero and we have 
the free Schrodinger equation 

^ dxl ~' df ^ ' 

That is the effect of the 5-functions are confined to the boundary of TZ, e.g. on the 
hyperplanes xj = xj+i, and their effect can be formulated as a boundary condition on the 
hyperplanes (we use also the fact that "^s is a symmetric function): 



d d 



dxj+i dxj 



IXj^l—Xj 



Thus the problem is to solve (JTj) in 71° subject to the boundary conditions ([H]) and the 
initial condition in TZ 

N 

^s{x;0) = l[5{x,-y,) (9) 

i=i 

where t/j satisfy ([5]) and we've dropped the normalization constant A^! since it can be 
incorporated at the end. The solution is given below in f|T4|) . 

2 Bethe Ansatz 

We now explain how the ideas of Bethe Ansatz [21 [7] can be employed to solve the time- 
dependent problem without using the spectral theory of the operator H. Since we avoid 
the eigenvalue problem, there are no Bethe equations in this approach^ To set notation 
and to see the argument in its simplest form, we first solve N = 2. 

2.1 N = 2 

It is elementary to verify that 

A{ki, A;2)e^(*^i^i+'=2^2)e-^(^(^^)+^(^2))* dkidk2 




with 

e{k) = e 

"^The methods here are an adaptation of the methods of pT| . 
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solves ([7]) in Xi < X2. The insight from Bethe Ansatz is to add to this solution another 
solution with the integration variables permuted, 

so that the boundary condition ([8]) can be applied pointwise to the above integrand. The 
result is the condition 

ik2Ai2 + ikiA2i - ikiAi2 - 'ifc2^2i = c{Ai2 + A21), 

or equivalently, 

c-i{k2-ki) 

M\ — --71 ^12- 

c + i[k2 - ki) 

We set 

.„.^s„.(*.-M--^^||^ (10) 

SO that the above reads A21 = 5*21^12. It is also convenient to set 

S(k) = -^^. (11) 

Observe that S{k) extends to a holomorphic function in the lower half-plane and so its 
Fourier transform, S{z), is supported in [0, 00). 

With the initial condition ([9]) in mind, we choose 

Ai2{ki, k2) = e''^^^y^+^^y^'^ 

so the above solution becomes 

\^^i{k,{x^-yi)+k2{x2-y2)) + S2l{k2, A;i)e*('=2(xi-y2)+fci(x2-yi))j ^-i{e^^+e^^)t ^J.^^).^ (^12) 

where from now on dkj dkj/27T. 

At t = the first term evaluates to 6{xi — yi)S{x2 — 1/2) ■ Thus we must show 

S{k2 - A;i)e*^i(^2-2'i)+*^^(^i-2^2) dkidk2 = 





for yi < y2 and Xi < X2 for the solution f[T2|) to satisfy the initial condition Qj. Making 
the change of variables ^2 ^ ^2 + ki and ki — ^ ki in the above integral, we have (after 
performing the resulting integration over ki) 

S{xi + X2-yi~ y2) I S{k2)e'^'^''^-y^^ dk2 = 5{xr +X2-yi- y2)S{x^ - 1/2). 



Now S{xi — ^2) is nonzero only in the region 

and the deha function requires that Xi + X2 = yi + 2/2 ior a nonzero contribution. But in 

2xi < xi + X2 = yi + y2 < 2?/2, 

i.e. xi < y2. Hence the integral in the region 71 is zero; and thus, we conclude (fT2|) is the 
sought after solution for N = 2. 

2.2 General N 

Let a E Sn he a permutation of {1, . . . , A^}. Recall that an inversion in a permutation a 
is an ordered pair {a{i),a{j)} in which i < j and cr{i) > cr(j). We set 

Aa- = Y\ {Sa/3 '■ {«5/3} is an inversion in a} (13) 

where Sa/s is defined by ([10]). Thus, for example, A231 = S21S31 and Aid = 1- We claim 
the solution that satisfies ([7]) with boundary conditions ([H]) and initial condition is 

^s{x; t)= J2 / ■ ■ ■ / n e'^^'^w^"^-^'-'^)^ e-''^^ '^^^^ dk,--- dkN- (14) 

o-g5]v j=l 

First, it is clear that f|T^ satisfies ([7]) in 7^°. As demonstrated in [7], if Tjcr denotes a 
with the entries a{i) and cr(z + 1) interchanged, the boundary conditions will be satisfied 
provided that 

for all a. Let us see why this relation holds. Let a = <j{i), P = cr{i + 1), and suppose 
a > [3. Then {a, (3} is an inversion for a but not for Tjcr, so S'q,/3 is a factor in Ag- but not 
in ^Ticr, and all other factors are the same. Therefore, using Sai^Spa = 1, we have 

The same identity holds immediately if (3 > a, since {(3, a} is an inversion for TjO" but 
not for a. 

As in the N = 2 case, the term corresponding to the identity permutation in 
satisfies the initial condition . Thus to complete the proof we must show 

N 



E 



/"■■■/ IT e^'^-w^^^'^'^'^f^)) dki--- dkN = 0. 
Jr Jr 



Let I (a) denote the integral corresponding to a term in the above sum. Recalhng the 
definition (fTTjl of S, the integrand for J(cr) becomes 

N 

J^e»fcaO)(^j-3/<.0)) JJ" [S{ka - kp) : {a, {3} is an inversion in a) . 

We again use the fact that S{k) extends analytically into the lower half-plane. 

If a number 7 appears on the left side of an inversion and never on the right side then 
in the integrand it appears only in factors of the form S{k^ — kp) and so integrating with 
respect to k^ shows that the integrand is zero unless 

a^a-i{7) > y-y- (15) 

Similarly, if a number 5 appears on the right side of an inversion and never on the left side 
then in the integrand it appears only in factors of the form S{ka — ks) and so integrating 
with respect to ks shows that the resulting integrand is zero unless 

Xa~^5) < ys- (16) 

We show I{a) = by induction on A^. If o"(l) = 1 then 1 does not appear in an 
inversion, we can integrate with respect to fci, and we are reduced to the case — 1. So 
assume a{l) = 7 > 1. Then (7, 1) is an inversion and there is no inversion of the form 
(a, 7) because 7 appears in slot 1. So we can apply (fT5l) with this 7 and get the resulting 
integrand is zero unless 

xi > y.,. 

Next we observe as before that (7, 1) is an inversion, but now that 1, which appears on 
the right side of the inversion, cannot appear on the left side of an inversion (obviously). 
So we can apply (ITB]) with 5 = 1 and get 

Since yi < y^ the two inequalities give Xo-i(i) < Xi, which cannot happen. This completes 
the proof of ffTI|) . 

The limit of impenetrable bosons is the limit c — > 00. In this limit Sa/s — > — 1 and 
ly^C'^) where 'mv{a) is the number of inversions in a. Thus in TZ 

*5(x; t) ^oo(x; t)= /"■■■/" det (e^'=-(^^-?'")-^*^'=-) dfci ■ ■ ■ dk^ (17) 

Jr Jr 

and the extension to follows by requiring \I'oo to be a symmetric function of Xj. 
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